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In case of the nonlinear Schrodinger equation with designed group velocity dispersion, variable 
nonlinearity and gain/loss; we analytically demonstrate the phenomenon of chirp reversal crucial for 
pulse reproduction. Two different scenarios are exhibited, where the pulses experience identical dis- 
persion profiles, but show entirely different propagation behavior. Exact expressions for dynamical 
quasi-solitons and soliton bound-states relevant for fiber communication are also exhibited. 

PACS numbers: 42.81. Dp, 05.45.Yv, 42.65.Tg 



Nonlinear Schrodinger equation (NLSE) is known to 
govern the pulse dynamics in nonlinear optical fibers [l| . 
In recent years, the study of nonlinear fiber optics, deal- 
ing with optical solitons, has attracted considerable at- 
tention since it has an important role in the development 
of several technologies of the 21st century Q. NLSE with 
distributed coefficients such as, group velocity dispersion 
(GVD), distributed nonlinearity and gain/loss, is being 
studied extensively in order to determine the effect of 
various distributed parameters on the pulse profile. 

In the realistic situation in a fiber, there arises non- 
uniformity due to variation in the lattice parameters of 
the fiber medium, as a result of which the distance be- 
tween two neighboring atoms is not constant through- 
out the fiber. It may also arise due to the variation of 
the fiber geometry e.g., diameter fluctuation. These non- 
uniformities influence effects such as, loss (or gain), phase 
modulation, etc, which can be modeled by making corre- 
sponding parameters space dependent. From a practical 
point of view, tailoring of various fiber parameters may 
lead to effective control of the pulse. This has been one of 
the prime motivation of a number of authors to analyze 
NLSE in a distributed scenario. 

Dispersion management (DM) has emerged as an im- 
portant technology to control and manipulate the light 
pulses in optical fibers 0, Q. Pulse compression has 
been demonstrated with appropriately designed GVD 
and nonlinearity in the presence of chirping 0, H, Q , as 
also through soliton effects Q. Adiabatic soliton com- 
pression, through the decrease of dispersion along the 
length of the fiber has been shown to provide good pulse 
quality The possibility of amplification of soliton 
pulses using a rapidly increasing distributed amplifica- 
tion with scale lengths comparable to the characteristic 
dispersion length has been reported @. It has been nu- 
merically shown that, in the case where the gain due to 
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the nonlinearity and the linear dispersion balance each 
other, equilibrium solitons are formed [Toj]. Serkin and 
Hasegawa have formulated the effect of varying disper- 
sion and other parameters on the soliton dynamics and 
have explained the concept of amplification of soliton 

& 

The formal structure of the Lax pair for the deformed 
NLSE has been studied dUEl. In a significant result, it 
has been numerically shown that in an appropriately de- 
signed dispersion profile chirped pulses can be retrieved 
through chirp reversal at a calculated location in the fiber 
[3] • The advantage of pre-chirping of the input pulse in 
overcoming soliton interaction and dispersive-wave gen- 
eration has been noted earlier [TEj . 

In the present Letter, we demonstrate analytically the 
phenomenon of chirp reversal of quasi-solitons with a 
designed dispersion profile. Very interestingly we find 
two possibilities of chirp reversal for which the disper- 
sion profiles are identical. However, they exhibit entirely 
different propagation behavior. In one case the motion 
is sinusoidal and in the other it shows pulse acceleration. 
The procedure to control pulse dynamics is also pointed 
out. Exact expressions for dynamical quasi-solitons and 
soliton bound-states relevant for fiber communication are 
exhibited. 

It is worth emphasizing that, exact solutions have 
played crucial role in demonstrating different pulse shap- 
ing techniques. The soliton solutions of NLSE or modifi- 
cations of the same has come in handy in studying these 
mechanisms. In the same light, finding exact solutions 
for general types of distributed scenarios will illustrate 
the subtle effects and interplay of various parameters on 
formation and propagation dynamics of light pulses. 

We develop a methodology to obtain self-similar soli- 
tary wave solutions of generalized NLSE model with 
varying nonlinearity, GVD, gain/loss and a confining 
oscillator which can be further modulated. One or 
few of these parameters can be switched off depend- 
ing on the situation at hand. It is shown that, this 
equation decouples into elliptic function equation and a 
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Schrodinger eigenvalue problem. This allows one to ana- 
lytically treat a variety of distributed scenarios, a few of 
which we explicate in the text. In the context of Bose- 
Einstein condensates the procedure to deal with variable 
coefficient NLSE in the absence of GVD has been carried 
out recently by the present authors [l6[ • GVD leads to a 
fundamentally new control parameter in the present case 
dealing with optical fibers. For example, a subtle inter- 
play of GVD and nonlinearity leads to a soliton bound- 
state as will be seen below. The effect of GVD, alternat- 
ing between normal and anomalous dispersion, on pulse 
profile is also discussed. 

For the purpose of analytic demonstration of chirp re- 
versal we start with a NLSE model with variable GVD, 
nonlinearity and loss/gain 

id x q(z, t) + ^ld tt q(z, t) + ~f(z)\q\ 2 q(z, t) + ig(z)q(z, t) = 

With q(z,t) = a(z)u(z,t) and a(z) = 
exp[— J* dz' g(z')\, one obtains, 



idM*',t) + ^Y^d tt u(z',t) +7(*')MV*',*) = o, ( 2 ) 
where z' = f* dz"a 2 (z"). 

Keeping in mind, pre-chirping and self-similar nature 
of the pulse we make use of the following ansatz, 



where z" = J* p(s)ds. 

The above equation has been numerically investigated, 
where a chirp-reversal was observed for quasi-soliton hav- 
ing a profile intermediate to a Gaussian and the funda- 
mental NLSE soliton These are stationary solutions 
obeying NLSE with an additional oscillator term, which 
explains the above profile. The exact solutions of Eq. 
(0) can be obtained, following the formalism developed 
inRef. 0: 



u(z', t) = VK^) v [p(z% A exp [iC(z')t 2 ] , (3) 

where, p and C are real functions of z' . 

Defining r = p(z')t, for preserving space-time identity 
one obtains 



dz' ^ 



where 



8r 



d e p 2 d 2 v(z f ,r) 
~2 dr~ 2 
K x r 2 iK^v 



-^z')p\v\ 2 v(z',r) 
~- 0, (4) 
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p z > + Cd e p 
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C z > + C 2 d e 
P 3 



(5) 



We now tailor the dispersion profile with Kq = and 
K\ = const: 



. dv d e p 2 d 2 v(z',r) 2 , K x r 2 

'dz 1 + ~2 dr*~ + 7 ^ )P,Zy| ^ ' T) = ~2~"' (6) 

In order to map the above equation to one with con- 
stant anomalous dispersion we assume the constraint 
d e p = 1 to obtain, 



.dv \d 2 v(z",r) . , AI ,o , / x Kit 2 



v(z", t) = ^A(z^)F[A(z"){t - A(z")}]e^"> T \ (8) 
where 

*(*", r) = ai (z") + b{z")T - \c{z")t 2 . (9) 

o-(i) 

Here, ai (z") = a lo - ^ /* .4 2 ( 2 )<fe, A(z") = 
J Z v{z)dz, which satisfies the following equation: 

dA 



— +c(z")A(z") = b(z"), 



(10) 



with the general solution, 



A(z") = { e M*")dz"y-i x 



dz" {b{z")e^ c ^ dz "^ 



(11) 



The parameter c(z") obeys Riccati equation: 

c z „-c 2 {z")=K 1 , (12) 

which can be exactly mapped to linear Schrodinger eigen- 
value problem. We also find the following consistency 
conditions: 



j(z') = j A(z')/Ao, b(z") = A{z") 
A(z") = A exp { /*" c(z)dz} , A > 0. 



(13) 



F obeys elliptic function equation in new variable T = 
A(z"){t- A(z")}: 



F"(T) - \F(T) + 2kF 3 (T) = 0, where k = (14) 

The twelve Jacobian elliptic functions satisfy above 
equation. These functions interpolate between the 
trigonometric and hyperbolic functions in the limiting 
cases [17]. Bright soliton solutions of the type v(z" ,t) = 
y^Aiz 77 ) sech[T/To]e^*"' r ) exist for k > 0, where T 2 = 
— Ao/jo and A = 1/2Tq , similarly kink- type dark solitons 
exist for 70 > 0. We further note that, with normal dis- 
persion one can obtain dark solitons for 70 < 0. It needs 
to be emphasized that, in the present approach the oscil- 
lator term leads to a dynamical chirp and modulates the 
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FIG. 1: The figure depicting variation of chirp parameters 
with z. The red curve shows the launching chirp and the 
blue curve shows the same in combination with the dynamical 
chirp. 



pulse profile. However, the pulse retains its fundamental 
NLSE soliton character in the scaled variable z" . 

Below, we examine the formation of bright quasi- 
soliton like excitations, exhibiting chirp reversal phe- 
nomenon. This is accomplished by appropriate tailoring 
of GVD and pre-chirping of the launching pulse. Com- 
bining Eq. ([5]) and constraint d e p = 1, yield the following 
expressions for the tailored dispersion and the chirping 
parameter, 



d! d 2 d 

c = f: and de d^ 



dz' 



(15) 



It is interesting to notice that the choice of the constant 
Ki, gives rise to two scenarios, having identical disper- 
sion and chirp profiles, but possessing entirely different 
pulse velocities. 

We list below some explicit examples, depicting a va- 
riety of novel control mechanism for pulse manipulation. 
Soliton pulses exhibiting chirp -reversal- Inspired from 
the numerical investigations of Kumar and Hasegawa 
[HI], we first consider in Eq. (fT5|) . K\ > 0, which is 
equivalent to a regular oscillator potential in Eq. Q. 
Dispersion profile in this case reads, 




FIG. 2: Pulse propagation having sinusoidal velocity profile. 
The plot shows broadening after chirp-reversal. 




d e (z') = cosh(5z')- 



-^sinh((5z / ),with(5= [Ki+C 2 (0)] 

(16) 

From above dispersion profile we compute launched 
chirp parameter C(z') and plot it together with the dy- 
namical chirp in Fig. [U From Eqs. ([5|) and ([T6|) 
it is clear that launched chirp profile changes sign at 
z(z' c ) = (1/8) tanh-^-C^)/*]. 

The expression for traveling quasi-soliton, propagating 
with velocity v(z") = Aq cos (Kqz") , reads 



FIG. 3: Soliton pulse acceleration with parameter K\ < 
implying an expulsive oscillator scenario. 

reversal, we notice that the pulse seems to broaden, as 
is clearly seen in Fig. [2l Hence, the pulse needs to be 
retrieved at this point. With the help of a normal dis- 
persive element such as a grating the original pulse can 
be recovered [l8| . 

In contrast to the above example, if we consider K\ < 
0, we obtain 5 = [i^-C^O)] 1 / 2 , which again leads to the 
same dispersion profile as that of the previous example, 
-,y 2 but with the pulse velocity v(z") = Aq cosh (Kqz") . The 
expression for the soliton profile in this case is, 



q(z,t) = a(z)\/p(z') sec(Kiz") sech [sec(i^iz // ){r 
x exp{z [C(z')t 2 + <S>(z",r)}}. 



q(z,t) = a(z) v / p(z f )sech(K 1 z ff )sech [&xAi(K 1 z"){t - A(z")}} 
x exp{z [C(z')t 2 + . (18) 

It is interesting to observe that, compared with the 
previous case pulse broadening is significantly reduced 
for the same parameter values. 

Soliton bound- states. - Starting from Eq. (pQ) without tai- 
loring the dispersion profile, we proceed to obtain self- 



) ^similar solutions assuming the ansatz solution of the type: 
(17) 



The presence of the dynamical chirp shifts the chirp- q(z,i) 
reversal location slightly away from z c . Just after chirp 



yfMzjF [A(z) {t - T(z)}} exp {G(z) + i$>(z, *)} . 

(19) 
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FIG. 4: Soliton bound-state in a medium with two dispersion 
regimes. 




z 



FIG. 5: Multiple soliton bound- state with oscillatory 
gain/ loss. 

The parameters appearing in the Eq. ([T9]) can be 
straightforwardly evaluated from the Eqns. ([TT|h ([T2]h 
(Q3J) and soliton profile can be obtained from Eq. (p~4j) . 

Below, we explicate some examples of spatial bound- 
states of solitons, arising from interplay of GVD, non- 
linearity and gain/loss. Below, Fig. [H depicts a two- 
soliton bound-state. This arises in a medium, where both 
anomalous and normal dispersion regimes are smoothly 
connected. In the presence of periodic gain/loss one ob- 
serves modulation in the bound-state profile as is shown 
in Fig. [5j 

In conclusion, we have obtained exact soliton solutions 
exhibiting chirp-reversal, while retaining their original 
profile, crucial for pulse recovery in fiber optics. Two 
different soliton sectors differing in propagation behavior, 
but with identical dispersion profiles, are analytically ex- 
hibited. We have outlined a general formalism for obtain- 
ing self-similar solutions of nonlinear Schrodinger equa- 
tion, in the presence of distributed coefficients, from 
which earlier scenarios follow as special cases [19]. It is 
shown that, this nonlinear system, involving pulse prop- 
agation with group velocity dispersion (GVD), variable 
nonlinearity, variable gain, exactly decouples into ellip- 
tic function equation and a Schrodinger eigenvalue prob- 



lem. This opens up a number of possibilities to take 
into account a wide class of distributed scenarios, in close 
conformity with the experimentally achievable situations. 
This incorporates a number of special cases dealt ear- 
lier, in the context of pulse compression. We find that, 
apart from compression, one can achieve control over the 
pulse velocity, pulse profile, through interplay of group 
velocity dispersion, nonlinearity, gain/loss. Formation of 
soliton bound-states in a medium with GVD alternating 
between normal and anomalous dispersion, is discussed. 
In the presence of oscillatory gain/loss profile, we find 
multiple bound state structure. 

We acknowledge many useful discussions with Prof. G. 
S. Agarwal. 
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